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Abstract. We investigate the evolution of the galaxy luminosity function from the VIMOS-VLT Deep Survey 
(VVDS) from the present to z = 2 in five ([/, B, V, R and /) rest-frame band-passes. We use the first epoch 
, - VVDS deep sample of 11,034 spectra selected at 17.5 < Iab < 24.0, on which we apply the Algorithm for 

Itn] ' Luminosity Function (ALF), described in this paper. We observe a substantial evolution with redshift of the 

global luminosity functions in all bands. From z = 0.05 to ^ = 2, we measure a brightening of the characteristic 
magnitude M* included in the magnitude range 1.8 — 2.5, 1.7 — 2.4, 1.2 - 1.9, 1.1 - 1.8 and 1.0 - 1.6 in the U, B, 
V, R and / rest-frame bands, respectively. We confirm this differential evolution of the luminosity function with 
rest-frame wavelength from the measurement of the comoving density of bright galaxies (Af < M*{z — 0.1)). 
This density increases by a factor of around 2.6, 2.2, 1.8, 1.5, 1.5 between z = 0.05 and z; = 1 in the U, B, V, R, 
I bands, respectively. We also measure a possible steepening of the faint-end slope of the luminosity functions, 
with Aa ~ — 0.3 between z = 0.05 and z = 1, similar in all bands. 
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1. Introduction 

The luminosity function (LF) of field galaxies is a funda- 
mental diagnostic of the physical processes that act in the 
formation and evolution of galaxies. The LF evolution is 
mainly determined by the combination of the star forma- 
tion history in each galaxy and the gravitational growth of 
structures, through merging. These two different processes 
are better probed by the luminosity emitted in the blue 
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and red rest-frame wavelengths, respectively. The relative 
contribution of these processes to the cosmic history is re- 
flected in the LF evolution, which therefore is expected to 
be different as a function of rest-frame wavelength. Large 
deep redshift surveys, combined with multi-color imaging, 
are necessary to perform this measurement. 

The local LF is now well constrained by the results 
of two large spectroscopic surveys: the Two-Degree Field 
Redshift Survey (2dFGRS; [Norberg et 31^002) an d the 
Sloan Digital Sky Survey (SDSS: IBlanton et al. E 003V 
These measurements of the local LF represent the lo- 
cal benchmark for all studies of the LF evolution. Up 
to 2: ~ 1 the Canada-France Redshift Survey (CFRS; 

et al. | L995) represents a sample of 591 spectroscopic 
redshifts of galaxies, from which it was demonstrated that 
the global LF evolves with cosmic time. [Lilly et al. i ;i995) 
showed that the evolution of the LF depends on the galaxy 
population studied. The LF of the red population shows 
few changes over the redshift range 0.05 < z < 1, while the 
LF of the blue population brightens by about one mag- 
nitude over the same redshift interval. Up to z ~ 0.6, 
the Canadian Network for Observational Cosmology Field 
Galaxy Redshift Survey (C N0C2; Lin ct al. 1999|l and the 
ESO-Sculptor Survey (ESS; |de Lapparent et al. 2003| ) de- 
rived the LFs per spectral type with spectroscopic red- 
shift samples of ~ 2000 and 617 galaxies, respectively. 
They confirmed a steep faint-end slope of the LF for the 
blue galaxy types. At higher redshift, LF measurements 
based on photometric redshifts have been derived by, e.g., 
IWolf et al. r 2003) up to z < 1.2. ICabasch et al. r 2004) up 
to z < 5. Samples of Lyman-break selected galaxies have 
also been used to measure the LF at such high redshift 
3 < z < 5 fe.g..lSteidel et'^999V 

The VIMOS (Visible Multi-Object Spectrograph) 
VLT (Very Large Telescope) Deep Survey (VVDS) 
is a deep spectroscopic survey conducted over a 
large area associated with multi-color photometric data 
l|Le Fevre et al. 2004all . Because of its characteristics, the 
VVDS is very well suited for detailed studies of the LF 
evolution : 

— the Universe is surveyed over more than 90% of its cur- 
rent age with spectroscopic redshifts, which allows us 
to measure the LF evolution in a coherent way within 
a single sample from z — 0.05 up to high redshift; 

— the spectroscopic targets are selected on the basis of a 
simple magnitude limit criterion, with no attempt to 
exclude stars or AGNs, or to select objects on the basis 
of their colors or morphology. Therefore, this selection 
minimizes any bias in sampling the galaxy population 
up to high redshift; 

— the multi-color coverage allows us to span a large range 
of rest-frame wavelengths, related to different physical 
processes, and thus to derive the LFs in several rest- 
frame bands; 

— the spectra can also be used to derive the LFs as a 
function of specific spectral properties. 



In this paper, we focus on the deep fields of the VVDS, 
which are the VVDS-0226-04 field IILe Fevre et al. £ 005) 
and the VVDS -Chandra Deep Field South (CDFS; 
ILe Fevre et"aT~b 004b') . The first epoch VVDS deep sam- 
ple contains 11,034 spectra of objects selected at 17.5 < 
Iab < 24.0. The first goal of the LF analysis is to charac- 
terize the statistical properties of the whole population. It 
is the first measurement that a theory of galaxy formation 
must account for, because it is free of any of the possible 
biases which can arise when the sample is split into dif- 
ferent sub-samples according to various selection criteria 
(e.g., rest- frame colors, morphology, spectral properties). 
In this paper, we present the measurement of the evolution 
of the global LF up to z = 2. To investigate the depen- 
dence of this evolution on the rest-frame wavelength, the 
global LFs are estimated in five rest-frame bands, which 
span the wavelength range 3600A< Ae// < 7840A. 

The paper is organized as follows. In Section 2 we 
briefly present the VVDS Deep first epoch sample. In 
Section 3 we describe the target sampling rate and the 
spectroscopic success rate of our data. In Section 4 we 
discuss two points relevant to the estimate of the global 
LF with VVDS data. This estimate is performed with our 
LF tool named Algorithm for Luminosity Function (ALF), 
extensively described in the Appendix. In Section 5 we 
present our results, compared with other literature mea- 
surements in Section 6. Conclusions are presented in 
Section 7. This paper will be followed by an analysis of the 
evolution of the LF per spectral type IjZucca et al. £ 005). 
and as a function of environment l|Ilbert et al . 2005). 

We use a flat lambda {ilm = 0.3, = 0.7) cosmol- 
ogy with h — i/o/100 km s^^ Mpc^^. Magnitudes are 
given in the AB system. 

2. Data description 

We consider the deep spectroscopic sample of the first 
epoch data in the VVDS-0226-04 and VVDS-CDFS fields. 

McCracken et al. (2003) describe in detail the pho- 
tometry and astrometry of the VVDS-0226-04 field ac- 
quired with the wide-field 12K mosaic camera at the 
Canada- France-Hawaii Telescope (CFHT). The deep field 
covers 1.2 deg^ and reaches the limiting magnitudes of 
Bab 26.5, Vab 26.2, Rab - 25.9 and Iab - 25.0, 
corresponding to 50% completeness. These data are com- 
plete and free of surface brightness selection effects at 
Iab < 24.0, corresponding to the limit of the VVDS spec- 
troscopic sample. Apparent magnitudes are measured us- 
ing Kron-like elliptical aperture magnitudes ( Kron 1980|l . 
with a minimum Kron radius of 1.2 arcsec. They are cor- 
rected for the galactic extinction estimated at the center 
of the VVDS-0226-04 field. For a large fraction of the field 
we have also U band data, taken at the ESO 2.2m tele- 
scope and reaching a limiting magnitude of Uab ^ 25.4 
(^Radovich et al. 2004|l . 

For the VVDS-CDFS, w e have used the EIS /-band 
photometry and astrometry (' Arnouts et al. 2001|l for our 
target selection, and the multi-color U, B, V, R, and 
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/ photometric catalogue from the COMBO-17 survey 
II Wolf et al. 200411 . 

The VVDS redshift survey uses the high multi- 
plex capabilities of the VIMOS instrument installed at 
the Nasmyth platform of Melipal of the VLT-ESO in 
Paranal (Chile). The spectroscopic observations were ob- 
tained during two runs between October and December 
2002. The spectroscopic targets were selected from the 
photometric catalogues using the VLT-VIMOS Mask 
Preparation Software (VVMPS: iBottini et al.'b OOS). The 
spectroscopic multi-object exposures were reduced us- 
ing the VIPGI tool dScodeggio et al. 2005| ). The sam- 
ple of spectroscopic redshifts obtained in the VVDS- 
CDFS is described in ILe Fevre erarT 2004b) and the 
sample obtained in the VVDS-F02 is described in 
ILe Fevre et al. r 2005^. A total of 11,034 spectra were ac- 
quired as primary targets in the two fields. The range of 
magnitude of the observed objects is 17.5 < Iab < 24.0. 
The deep spectroscopic sample (VVDS-0226-04-HVVDS- 
CDFS) consists of 6582-^1258 galaxies, 623-^128 stars and 
62-1-9 QSOs with secure spectroscopic identification, i.e. 
quality flags 2, 3, 4 and 9 (flags 2, 3, 4 correspond to 
redshifts measured with a confidence level of 75%, 95%, 
100%, respectively; flag 9 indicates spectra with a single 
emission line). 1439-1-141 objects have an uncertain red- 
shift measurement, i.e. quality flag 1 (flag 1 corresponds 
to a confidence level of 50% in the measured redshift). 
690-1-102 objects have no spectroscopic identification, i.e. 
quality flag 0. This sample covers ^ 1750 + 450 arcmin^, 
with a median redshift of about 0.76. The la accuracy 
of the redshift measurements is estimated at ~0.001 from 
repeated VVDS observations l|Le Fevre et al. E 005V 

3. Treatment of unidentified sources 

In the estimate of the luminosity function, we introduce 
a statistical weight Wi, associated with each galaxy i with 
a secure redshift measurement. This weight corrects for 
the non-observed sources and those for which the spec- 
troscopic measurement failed (unidentified sources). This 
method yields the best statistical estimate of the total 
number of galaxies with the same properties as galaxy i, 
in the full field of view sampled by the spectroscopic data. 
The statistical weight Wi is the product of : 

— the weight wf^^, inverse of the Target Sampling Rate 
(hereafter TSR). The TSR is the fraction of objects 
in the photometric catalogue that have been spectro- 
scopically targeted. It can be a constant or a function 
of a number of parameters according to the strategy 
adopted in selecting the spectroscopic targets. 

— the weight wf^^, inverse of the Spectroscopic Success 
Rate (hereafter SSR). The SSR is the fraction of the 
spectroscopically targeted objects for which a secure 
spectroscopic identification has been obtained. In this 
paper, the LFs arc computed using galaxies with secure 
redshift measurements, i.e. with spectroscopic quality 
flag 2, 3, 4 and 9. The SSR is the ratio between the 



number of objects with high quality flag 2, 3, 4 and 
9 and the total number of spectra (quality flags 0, 1, 
2, 3, 4 and 9). The correct estimate of the SSR is not 
trivial, because it can be a function of a large number 
of parameters, like, for example, magnitude, surface 
brightness, redshift and spectral type of the objects. 

3.1. The TSR and its associated weiglit w'^^^ 

The VVDS strategy in selecting spectroscopic targets has 
been to select targets quasi-randomly from the photo- 
metric catalog, thus minimizing any bias in sampling the 
galaxy population. In the random selection process, the 
VVMPS tool HBottini et al. B 005) uses the information 
about the size of the objects in order to maximize the 
number of slits per VIMOS pointing. As a consequence, 
the flnal spectroscopic sample presents a bias with respect 
to the photometric one, with large objects being under 
represented fsee lBottini et al"^ 005. for a discussion). The 
parameter used by VVMPS to maximize the number of 
slits is the x-radius, which is the projection of the angular 
size of each object on the x-axis of the image, correspond- 
ing to the direction in which the slits are placed. The 
x-radius is deflned as x-radius=(n -I- 0.5) x 0.204, where 
0.204 is the pixel size of the image expressed in arcsec 
and n is an integer corresponding to the size of the ob- 
ject in pixels. The TSR in the VVDS-0226-04 is shown 
as a function of the x-radius in the top panel of Fig^ 
The TSR runs from ~ 25% for the smallest objects, to 
^ 10% for the largest ones. As shown in the bottom panel 
of Fig^ a large fraction of the total population (~ 75%) 
is targeted with TSR ~ 25%. The under-sampling of the 
largest objects (x-radius > 1.7) concerns less than 4% of 
the total population, targeted with TSR < 15%. Since the 
x-radius is the only parameter used to maximize the num- 
ber of slits, the correction to be applied in order to correct 
for this bias is well deflned and corresponds to using the 
weight wf^^ — l/TSR{ri), where is the x-radius of the 
galaxy i. 

3.2. The SSR and its associated weight w^^^ 

The second weight to be used in the estimate of the LF 
is wf^^, which is the inverse of the SSR. In Fig. 16 of 
ILe Fevre et"aO 2005l it is shown that the SSR is, as ex- 
pected, a function of the Iab apparent magnitude. The 
SSR is greater than 90% for Iab < 22.0 and smoothly 
decreases down to ^ 70% in the faintest half a magni- 
tude bin. In a flrst approximation, we could use this SSR 
distribution to derive wf^^^ as a function of the Iab ap- 
parent magnitude. However, this procedure implies that 
the objects with quality flags and 1 belong to the same 
population of the objects with a secure spectroscopic iden- 
tiflcation (flags 2, 3 and 4). In particular, it implies that 
they have the same redshift distribution. The redshift dis- 
tributions of galaxies with quality flags 4, 3 and 2 are 
shown in Fig [3 The distributions for each flag are clearly 
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Fig. 1. Top panel: Target Sampling Rate as a function of 
the x-radius for the VVDS-0226-04 field. Bottom panel: 
x-radius distribution in the VVDS-0226-04 photometric 
catalogue at 17.50 < Iab < 24.0. 



different, reflecting the fact that the quality flag is related 
not only to the signal-to-noise of the spectrum, but also to 
the number and the strength of identifiable spectral fea- 
tures. This suggests that the galaxies with quality flags 
and 1 are likely to have a different redshift distribution. If 
that is the case, we would not be allowed to use wf^^ as 
a function of magnitude only. 

Therefore, making use of the multi-color properties of 
our sample, we have analyzed the distribution of photo- 
metric redshifts for the spectroscopic targets with flag 
and 1. For this analysis we have used only a subset area 
of the VVDS-0226-04 field, with - 1100 spectra, in which , 
in addition to the U photometry ( Radovich ct al. 2004), 
we have also J and K photometry (l ovino et al. 2005 ). 
We have restricted this analysis to the area with near- 
infrared data, since near infrared photometry allows us 
to estimate robust photometric redshifts at least up to 
z ^ 2 (see IBolzonella et al. 20051 for a detailed descrip- 
tion of the method). A redshift probability distribution 
function (hereafter PDFz) is estimated for each object of 
the spectroscopic sample, using the photometric redshift 
code of Le Phare^ (Arnouts & Ilbert). We sum the normal- 
ized PDFz of all galaxies to estimate the expected redshift 
distribution l|Arnouts et al. 2002|l . The stars are removed 
from the sample on the basis of their spectral identifica- 
tion if they have high quality spectroscopic flags, or on 
the basis of photometric criteria for the low quality flags 
(jBolzonella et al. 2005(1 . The estimated redshift distribu- 
tion of galaxies with quality flag and 1 is shown in the 
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Fig. 2. Redshift distributions normalized to unity for each 
spectroscopic quality flag. For high quality flags 2, 3, 4, we 
show both the spectroscopic redshift distributions (dashed 
lines) and the photometric redshift distributions (solid 
lines) on the same area. For quality flags and 1, we 
show only the redshift distribution estimated using the 
photometric redshifts. 



bottom right panel of Fig[21 As expected, the estimated 
redshift distribution of the low quality flag galaxies differs 
from the redshift distribution of high quality flag galax- 
ies, while the distributions of the photometric and spec- 
troscopic redshifts for galaxies with flag > 2 are consistent 
with each other. 

We have then derived the SSR in various bins of appar- 
ent magnitude as the ratio between the estimated redshift 
distribution of high quality flag galaxies (quality flags 2, 3, 
4, 9) and the estimated redshift distribution of all galax- 
ies (quality flags 0, 1, 2, 3, 4, 9). This SSR is shown in 
Fig'Elas a function of redshift in four apparent magnitude 
bins for Iab > 22.5 (at brighter magnitudes, the SSR is 
close to unity). Fig|31 clearly shows that the global SSR 
indeed decreases for fainter apparent magnitude bins and 
it varies significantly with redshift. The shape of the SSR 
is similar in all magnitude bins showing a maximum ef- 
ficiency in the redshift measurement at z ~ 0.7 and a 
minimum SSR for z < 0.5 and z > 1.5. The dependence 
of the SSR on the redshift is related to the presence of the 
[O II] line, and/or the Balmer break within the observed 
spectral window 5500A< A < 9500A. The weight w^^^ is 
shown in the bottom of each panel in Fig|21 The weight is 
binned in redshift in order to limit the statistical noise. 

At z > 1.5, the uncertainties on our weight are large 
due to the smaller number of galaxies, and to the un- 
certainties on the photometric redshifts at such redshifts 
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Fig. 3. Spectroscopic Success Rate as a function of red- 
shift and per apparent magnitude bin. The associated 
weight is shown in the bottom of each panel. 



IjBolzonella et al. 2005|l . We could also perform an other 
estimate of the weight at such high redshifts, using the 
spectroscopic redshifts of the spectra with quality flag 1 
(50% of confidence level) and assuming that quality flag 
objects have the same redshift distribution. Applying 
this method, we find SSR^20-30% above z > 1.5, which 
provides a weight ~2-3 times greater than the weight esti- 
mated with the photometric redshift method. The value of 
y^ssR at z > 1.5 will be refined in future analysis (Paltani 
et al., in preparation), using simulations and new spectro- 
scopic VIMOS observations with a blue grism. 

We apply the weights derived from this analysis to 
the whole sample, making the assumption that the subset 
area, from which they have been derived, is representative 
of our ability to measure a redshift. We assign to each 
galaxy, a weight wf^^ that depends both on the apparent 
magnitude rrii and on the redshift Zi of the considered 
galaxy i. 

To summarize, we have derived our statistical weights 
as the product of Wi — wf^^ x wf^^ = l/TSR{ri) x 
l/SSR{mi, Zi). This weighting scheme allows us to correct 
for: 

— the TSR, taking into account the small under-sampling 
of the largest objects in the target selection; 

— the SSR, taking into account the dependence on the 
apparent magnitude and redshift. 



4. ALF applied to the VVDS 

We have measured the LF on the VVDS data, using 
our luminosity function VVDS tool, named Algorithm for 



Luminosity Function (ALF). The methods implemented 
in ALF are extensively described in the Appendix. In this 
Section we briefly discuss two points which are relevant 
for a better understanding of our treatment of the data in 
this paper. 

4.1. Effect of the w^^^ weight on the LF estimates 

As shown in Section 3, while the weight w'^^^ is fully un- 
derstood and well established, the derivation of the weight 
^SSR jg jggg direct and subject to more uncertainties. To 
quantify the effect of w^^^ on our LF estimate, we have 
also derived an 'unweighted' LF, in which no correction for 
the SSR is appHed (i.e. w^^^ = 1). The 'weighted' and 
'unweighted' LFs arc shown in Fig^ in the B rest-frame 
band. 

Since the galaxies with flag and 1 can not be ignored 
in the LF estimate, the 'unweighted' LF is by deflnition a 
lower limit of the LF. However, given the relatively small 
fraction of galaxies with flag and 1, the difference in 
the overall normalization of the two LFs is small. From 
FigEl we see that the main effect of w^^^ is to steepen 
the slope of the 'unweighted' LF. In all the redshift bins, 
the parameter a of the weighted LF is smaller (i.e. steeper 
slope), by less than 0.2 up to z = 1.0, by less than 0.3 in 
the two higher redshift bins. This effect is clearly expected, 
since the galaxies with flag and 1, which are included 
in the LF estimate through the weight, are mainly faint 
galaxies close to our magnitude limit. Since a and M* 
are correlated, the steepening of the slope with the weight 
produces also a brightening of M* , less than 0.25 up to 
z = 2. 

4.2. Estimate of the LF from a homogeneous galaxy 
population 

Ilbert et al.(2004) have shown that the estimate of the 
global LF can be biased, mainly at its faint-end, when the 
band in which the global LF is measured is far from the 
rest-frame band in which galaxies are selected. This is be- 
cause different galaxy types have different absolute magni- 
tude limits, because of different k-corrections. In each red- 
shift range, we avoid this bias in our estimates of the LF by 
using only galaxies within the absolute magnitude range 
where all the SEDs are potentially observable. We perform 
only the l/Vmax estimate on the whole absolute magni- 
tude range. This estimator leads to an under-estimate of 
the LF in the absolute magnitude range fainter than this 
'bias' limit Ijllbert et al. E 004'). providing a lower limit of 
the LF faint-end. 

5. Results 

The global LFs are computed up to z = 2 in the flve stan- 
dard bands U, B, V, R, I (U Bessel, B and V Johnson, R 
and / Cousins). The LFs are computed using the weight- 
ing scheme described in Section 3. The U-, B-, V-, R- and 
/-band LFs are displayed in Figl^and in FigEl In Figd 
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Fig. 4. Comparison between the 'unweighted' LFs (sohd triangles for 1/Vmax, solid lines for STY) and our best LF 
estimates ('weighted') (open circles for 1/Vmax, dotted hues for STY) in the B rest-frame band. The best STY fits 
and the associated 90% error contours are shown as insets with the same symbols. We report in each inset the SDSS 
local estimate with open triangles. 



we plot the STY best estimate of a — M*, and the as- 
sociated error contours. For each band and each rcdshift 
bin, the Schechter parameters and the corresponding one 
sigma errors measured with the STY estimator are listed 
in Tabled 

5.1. Local LF at z = 0.1 

The local LF derived with the VVDS sample refers 
to the redshift bin [0.05 — 0.2]. The average redshift 
in this bin (< z >~ 0.14) is directly comparable to 
the average rcdshift of galaxies in local surveys with a 
brighter limiting magnitude, like the SDSS (< z >~ 
0.1). Due to the bright apparent magnitude cutoff of 
the VVDS sample {Iab > 17.5), the M* parameter 
of the STY fit in this redshift bin is essentially uncon- 



strained. Therefore, we set the M* parameter to the lo- 
cal value derived bv lBlanton et al."^ 2003V The LFs of the 
SDSS are expressed in the bands °-^u, °'^g, ° ^-^z 
HFukugita et al. 1996| ) blue-shifted at z — 0.1, which cor- 
respond roughly to the U, B, V, R, I bands of our stan- 
dard system. In order to check if the absolute magnitudes 
estimated in the SDSS band system and in our standard 
band system are comparable, we have estimated the ab- 
solute magnitudes in the filters °-^(/, °-^r, "'^i and 
from the apparent magnitudes measured in the in- 
strumental system (using the formulae A.l and A. 2 given 
in appendix A). The average difference between the ab- 
solute magnitudes computed in B, V, R, I and in 
0.1^^ o.ij,^ 0.1^ bands are less than 0.05. The difference is 
more significant in the U band (AM 0.25). We have 
therefore converted Mgjjgg{°-^u) to our band with the re- 
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Fig. 5. Estimate of the global LF in the U band from z = 0.05 to z = 2. The estimate is derived using the weighting 
scheme described in Section 3. We adopt the following symbols for the various estimators: circles for the l/Vmax, 
triangles for the SWML, squares for the C+ and solid lines for the STY. The vertical short-dashed lines show the faint 
absolute magnitude limits considered in the STY estimate (see Section 4.2). In each panel, we show also the local LF 
derived by the SDSS (long dashed hues). The vertical dot-dashed line in the redshift bin 0.05-0.2, corresponds to the 
faint absolute magnitude limit surveyed by the SDSS. 



lation M^jjssiU) = M^jjss{°-^u) -0.25. The local values 
of a and 0*, with M* set to the SDSS value, are listed in 
Tableland the LFs are shown in FiglSjand Fig|Bl Even 
if the volume surveyed by the VVDS in the first redshift 
bin is approximately one thousand times smaller than the 
volume surveyed by the SDSS, the estimates of the lo- 
cal LFs produced by the VVDS and the SDSS are in good 
agreement in the magnitude range in common to both sur- 
veys. However, in all the bands, the VVDS best fit slope 
is steeper than the SDSS slope. The larger difference is in 
the B band, where it is formally significant at ^ 3ct level 
(asDSS = -0.89 ± 0.03 while avvDS = -1-09 ± 0.05). 
In this band the VVDS slope is instead consistent with 
that derived by |Norberg et al. [ 2002) from the 2dFGRS. 



Even if the number of objects in the VVDS is smaller 
than in the SDSS, the faint-end slope of the LF is bet- 
ter constrained by the VVDS because it samples the lo- 
cal galaxy population about 3-4 magnitudes deeper than 
the SDSS. The steeper slope observed in the VVDS can- 
not be due to the effect of the applied weights since also 
the 'unweighted' LF, which under-estimates the slope (see 
Section 4.1), has a steeper best fit slope than the SDSS 
{avvDS — —1.02 ± 0.05 in the i3-band). The inclusion 
of a fit for simple luminosity and number evolution in 
the LF estimate, using the maximum likelihood estimator 
proposed bv IBlanton et al."T 2003). could also produce a 
flatter slope. 
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Fig. 6. Same symbols as in Fig[Sl in the rest-frame band B (upper-left), V (upper-right), R (lower-left) and / (lower- 
right). 



5.2. LF evolution up to z = 2 

The VVDS allows us to quantify the galaxy evolution 
within a single sample and with the same selection func- 
tion, over a wide redshift range. From z — 0.05 up to 
z = 2, the evolution of the bright part of the LFs is clearly 
evident from all non-parametric estimators shown in Fig|Sl 
and Fig|Sl It also appears to be a function of the consid- 
ered rest-frame wavelength. This can be quantified using 
the Schechter parameters measured with the STY estima- 
tor, as done below. 



To quantify the strength of the evolution with redshift, 
we have derived the density of galaxies brighter than the 
corresponding local value of M* : 

/^'^SDS S 
-oo 

where ^sdss local value estimated by 

Blanton et al. r 2003). This measurement quantifies 
the global evolution of the bright part of the LFs, 
in shape and in normalization. In all bands and up 
to z < 1, Mgjjgg is brighter than the faintest limits 
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Fig. 7. M* — a error contours, at 90% confidence level {2AlnC = 4.61), obtained with the STY estimator. The solid 
squares are the STY best estimates. Error contours are not drawn when one of the two values has been fixed (see text). 
When a is fixed, we estimate the range of allowed M* values, varying a between a = —2.0 and a = — 1.2. When M* 
is fixed, we vary M* by 0.5 around the set value. The error contours are shown for tlie U, B, V, R and / band LFs 
from the left to the right panels, respectively. The panels from the top to the bottom correspond to the eight redshift 
bins, from z = 0.05 to z = 2. We report with open triangles the SDSS local values (the error bars are included inside 
the symbols). 



used for the STY estimate. We have therefore limited 
this analysis to z < 1, in order to avoid extrapola- 
tion of the LF beyond the last measured point. In 
Fig|Sl wc show the evolution with redshift of the ratio 
p{M < Mgjjgg)/psDSs{M < M g ^ g g) foT thc five bands. 
As the figure clearly shows, the density evolution of 
bright galaxies is significantly dependent on the rest- 
frame wavelength, being stronger at shorter wavelengths. 
In the U band, p{M < Mgj^gg) increases continuously 
and becomes ^ 2.6 times greater than locally at z = 1.0, 
while in the / band this factor is only ~ 1.5 at z = 1. 

The evolution of the best fit M* as a function of red- 
shift for the five bands is shown in the central panel of 



Fig|3 We find that the characteristic magnitude M* of 
the whole population strongly evolves. Up to z = 1, the 
slope can still be constrained reasonably well and wc mea- 
sure a brightening of 1.57 ± 0.26, 1.48 ± 0.17, 1.41 ± 0.22, 
1.49 ± 0.25 and 1.45 ± 0.26 magnitudes in the U, B, V, 
R and / rest- frame bands, respectively. Above z = 1, the 
slopes are set to the a value obtained in the redshift bin 
0.8 < z < 1 and we measure a brightening of about 2.0, 
1.8, 1.4, 1.3, 1.2 magnitudes up to z = 2. When a is fixed, 
we estimate the range of allowed M* values, varying a 
between two extreme values of the slope, a = —2.0 and 
a = —1.2. We find a brightening included in the range 
1.8 - 2.5, 1.7 - 2.4, 1.2 - 1.9, 1.1 - 1.8 and 1.0 - 1.6 
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Fig. 8. Evolution with rcdshift of the ratio between the 
density of galaxies brighter than Mgjjgg{z = 0.1) and the 
SDSS local value. 

magnitudes in the U, B, V, R and / rest-frame bands, 
respectively. Also in this representation, the evolution is 
stronger in the bluer rest-frame bands. Since M* and a 
are correlated and we have some evidence that also a 
is changing with rcdshift (see below and upper panel of 
FigEJi we have verified that the observed evolution in A/* 
is not induced by the change with redshift of the a value. 
The bottom panel of FigEl shows the best fit M* param- 
eters derived by setting the value of a to the VVDS local 
value over the entire redshift range. Also in this case, a 
significant and differential evolution of M* is seen, with 

AM* 1.7, -1.6, -1.2, -1.1, -1.0 up to z ~ 2, in the 

U, B,V, R, I bands respectively. The measurement of this 
brightening is slightly sensitive to the adopted weighting 
approach and it is also measured, at a similar level, with 
the 'unweighted' LFs (see Fig^. 

The upper panel of FiglHl shows the best fit values of 
a as a function of redshift. The one sigma error bars on 
a take into account the correlation between a and M*. 
The data suggest a steepening of the slope with increasing 
redshift. The measured variation of a between z — 0.05 
and 2; = 1 is Aa ^ —0.3, similar in all the bands. 

6. Comparison with previous redshift surveys 

6.1. Comparison with the CFRS survey 

Lilly et al. (1995) have derived the global S-band LFs of 
the Canada- France Redshift Survey (CFRS) up to z ~ 1. 
The CFRS spectroscopic sample contains 591 redshifts of 
galaxies selected with 17.5 < Iab < 22.5. The survey cov- 
ers 125 arcmin^ in five separated fields. The VVDS deep 




1+Z 

Fig. 9. Evolution in the five bands of the parameter a 
in the top panel and of the parameter M* in the second 
panel. The bottom panel shows the evolution of M* with 
the slope fixed to the VVDS local value. 



spectroscopic sample is surveying the galaxy population 
1.5 magnitude fainter and the field of view is ten times 
larger than the CFRS. The comparison between VVDS 
and CFRS results in three redshift bins is displayed in 
FigEl (in the cosmology Q,m = 1, f^A = 0, which was 
adopted in the CFRS analysis). The estimated LFs for the 
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MAB(B)-51og(h) 

Fig. 10. Comparison between the CFRS and the VVDS 
global B-band LFs. The solid lines (STY) and the circles 
(1/Vmax) correspond to the VVDS estimates. The vertical 
short-dashed lines are the faint absolute magnitude limits 
considered in the STY estimates. The long-dashed lines 
correspond to the CFRS estimate. The vertical dot-dashed 
lines correspond to the faint absolute magnitude limits 
surveyed by the CFRS. 

two surveys are in agreement up to the faintest absolute 
magnitude limits reached by the CFRS. The slopes of the 
VVDS are, however, steeper than the CFRS slopes (the 
difference is particularly significant in the redshift bin ]0.5- 
0.75]). The slopes estimated from the VVDS are clearly 
more robust, since the VVDS is 1.5 magnitudes deeper 
and contains 10 times more galaxies than the CFRS. 




M,B(B)-51og(h) 

Fig. 11. Comparison between the global B-band LFs 
derived with the HDF data UPoli et al. E 003) and with 
the VVDS data. The solid lines (STY) and the circles 
(1/Vmax) correspond to the VVDS estimates. The vertical 
short-dashed lines are the faint absolute magnitude lim- 
its considered in the STY estimates. The long dashed lines 
correspond to the global LFs derived bv lPoli et al. f 2003'). 
The vertical dot-dashed lines correspond to the faint ab- 
solute magnitude limits surveyed by the HDF data. 

6.2. Comparison with the HDF data 

Poli et al. (2003) have derived the global i?-band LFs from 
z = 0.4 up to z = 3.5, using a composite sample of 1541 
I-selected galaxies down to Iab — 27.2 and 138 K-selected 
galaxies down to Kab — 25. The faintest galaxies of this 
composite sample are drawn from HDF North and South 
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data. Data from two additional fields (the CDFS and the 
field around the QSO 005 5-269), on which th e K20 spec- 
troscopic survey is based ICimatti et al. 2002|l . have been 
added to constrain the bright-end of the LF. Given the 
faintness of this sample, most of the redshifts (~ 80%) are 
photometric redshifts. IPoli et arT 2003l have derived the 
global i?-band LF using the /-selected sample up to z = 1 
and the JsT-selected sample for 2; > 1.3. The HDF data sur- 
vey the LF faint-end about 2-3 magnitudes fainter than 
the VVDS data. The LFs from IPoh et"a rT2003'l and the 
corresponding VVDS LFs are shown in Fig^Jin three red- 
shift bins. At z > 1.3. IPoli et al. r 2003) have derived the 
LF in the redshift bin [1.3-2.5], that we compare here with 
our measurement in the redshift bin [1.3-2]. As shown in 
FigEl there is an excellent agreement in the bright part 
of the LF between the VVDS and the IPoh et al. r 2003) 
measurements, up to z = 2. In the faint part of the LF, 
the slope estimated bv lPoli et al. f 2003) is slightly steeper 
(Aa ~ 0.15) than the slope estimated with the VVDS 
data in the redshift bin [0.4-0.7]. 



6.3. Comparison with the COMBO-17 survey 

Wolf et al. (2003) have derived the LFs up to z = 1.2 with 
a sample of ~ 25,000 galaxies from the COMBO-17 survey. 
This sample is selected in the R band [Rvega < 24). The 
redshifts are photometric redshifts derived from medium- 
band photometry in 17 filters. The Schechter parameters 
of the COMBO-17 global LF are available in the online 
material of the paper llWolf et al. 003^1 . The comparison 
between the B-band global LFs of VVDS and COMBO-17 
surveys is shown in the Fig^|in five redshift bins up to 
z = 1.2. 

The bright parts of the LFs appear to be roughly in 
agreement, although some significant differences are seen 
in a few redshift bins (see, for example, the redshift bins 
[0.4-0.6] and [0.8-1]). Given the errors on the a - M* 
parameters reported by the two surveys, the overall LF 
shapes are not consistent with each other (see insets in 
FigEJ- Since the COMBO-17 sample is selected from the 
R band, its global i3-band LF could be affected by the bias 
described in Ilbert et al. f 2004) at z > 0.5. This bias intro- 
duces an overestimate of the LF faint-end at z > 0.5 and 
could explain the significantly steeper slope measured by 
the COMBO-17 survey in the redshift bins [0.6-0.8] and 
[0.8-1]. Since a and M* are correlated, the same effect 
could also explain the differences seen in the bright part 
of the LF. These discrepancies can also be due to other 
reasons as, for example, a smaller fraction of very blue 
galaxies in the I-selected VVDS sample (in fact, the LF of 
the bluest galaxies has the steepest faint-end LF) or a bias 
in the COMBO-17 estimate due to their use of photomet- 
ric redshifts. These possibilities will be better investigated 
through a comparison of the COMBO-17 and VVDS LFs 
for each galaxy type l|Zucca et al. E 005). since such a com- 
parison is much less affected by the bias discussed above 
(|Tlhert et al. E 004V 



7. Conclusions 

We use the first epoch spectroscopic deep sample of the 
VVDS, with 11,034 spectra selected up to Iab = 24, to 
derive the global LF up to z = 2 in the five bands U, B, 
V, R, I. The global LFs are measured using ALF and care 
is taken to remove the bias introduced by the difference 
of visibility of the different galaxy spectral types. 

We observe a clear evolution of the global LF with 
redshift in all bands and we find that this evolution is sig- 
nificantly dependent on the rest-frame wavelength, being 
stronger at shorter wavelengths. The comoving density of 
the bright galaxies increases with redshift from z = 0.05 
up to z = 1. This increase is by a factor 2.6 in the 
U band and becomes smaller for redder rest-frame wave- 
lengths, with values of the order of 2.6, 2.2, 1.8, 1.5, 1.5 in 
the U, B, V, R, I bands, respectively. 

In order to better distinguish the processes responsi- 
ble of this evolution, we have studied the evolution with 
redshift of the Schechter parameters computed with the 
STY estimator. This analysis suggests a possible steep- 
ening of the slope with redshift. The observed change 
in a is ^ —0.3 from z = 0.05 up to z = 1, similar in 
all bands. This evolution has to be confirmed with the 
on-going second epoch VVDS data, which will allow us 
to decrease significantly the statistical errors on a. This 
evolution of the global LF slope is expected because of 
the different evolutions observed for the different galaxy 
types l|Zucca et al. £ 005*). In particular, since the LF of 
blue galaxies has a steep slope and evolves strongly with 
redshift (e.g., [Lilly e'rar| l995, IZucca et al.'b 005;i. the rel- 
ative contribution of the blue population to the global LF 
increases with redshift and could explain the steepening 
of the slope. 

We also measure a significant brightening of the global 
LF with redshift. This brightening, parameterized as the 
change of the best fit value of M* , is a function of the rest- 
frame wavelength. Compared to the local SDSS values, 
we obtain a brightening included in the range 1.8 — 2.5, 
1.7 - 2.4, 1.2 - 1.9, 1.1 - 1.8 and 1.0 - 1.6 magnitudes 
from z = 0.05 up to z = 2, in the U, B, V, R and 
/ rest-frame bands. This tendency of a stronger bright- 
ening toward bluer rest-frame wavelengths is consistent 
with existing measurements at shorter and longer rest- 
frame wavelengths. In the rest-frame far-UV (1530 A), 
Arnouts et al. (2005), measure a brightening AM* ^ —2 
magnitudes up to z = 1, stronger than our measure- 
ment in the U band in the same redshift interval. In the 
near-IR, (Pozzetti et al. (2003) measure an evolution con- 



sistent with a mild luminosity evolution both in the J 
and K bands with AMj ~ -0.7 and AMk ~ -0.5 at 
z ~ 1. This differential evolution of M* with wavelength 
is expected, since the rest-frame luminosity at different 
wavelengths probes different physical processes acting in 
galaxy formation and evolution. The fact that the bright- 
ening is stronger in the bluest bands suggests that most 
of the evolution of the global LFs up to z = 2 is re- 
lated to the star formation history, better probed with 
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Fig. 12. Comparison between the COMBO-17 and the VVDS global B-band LFs. The solid lines and the points 
correspond to the VVDS estimates. The vertical short-dashed lines are the faint absolute magnitude limits considered 
in the STY estimates. The long dashed lines are the global LFs derived bv' Wolf et al. r 2003). The vertical dot-dashed 
lines correspond to the faint absolute magnitude limits surveyed by the COMBO-17 data. The best estimated values 
for the a, M* parameters measured by the VVDS are shown with solid squares in insets, and the associated 90% error 
contours with solid lines. The (a, M*) parameters given by COMBO-17 are the open stars, with error bars included 
inside the symbols. 



the luminosity measured at short rest-frame wavelengths. 
The luminosity density and star formation rate derived 
from the VVDS first epoch observations will be presented 
in ITresse et al. 1 20051 We will explore the evolution of 
the LFs per spectral types and as a function of the lo- 
cal environment in forthcoming papers IjZucca et al. E 005. 
Illbert et al. fe 005l 
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(a) Number of galaxies in the redshift bin (sample used for 1/Vmax estimate) 

(b) Number of galaxies brighter than bias limit (sample used for STY, C-f-, SWML estimate) 



Table 1. Schechter parameters and associated one sigma errors {2AlnC = 1) of the global LFs between z = 0.05 and 
2 = 2 and derived in the U, B, V, R, and / filters of the standard system. Parameters listed without errors are set 
'ad-hoc' to the given value. 
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Appendix A: The Algorithm for Luminosity 
Function (ALF) 

This section describes the standard methods implemented 
in our Algorithm for Luminosity Function (ALF) devel- 
oped within the VVDS framework. We present how we 
derive the rest-frame absolute magnitudes and the details 
of the 1/Vmax, C+, SWML and STY estimators imple- 
mented in this tool. 



A.l. Absolute magnitudes 

The k-correction depends on the galaxy spectral en- 
ergy distribution (SED). At high redshift, it is the main 
source of error and systematic in the absolute mag- 
nitude measurement. Using Le Phare, we adjust the 
best SED template on [/, B, V, R and I apparent 
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Fig. A.l. Difference between our 'rederived' absolute 
magnitudes and 'true' absolute magnitudes from GALICS 
as a function of redshift, in the U, B, V, R and / reference 
bands from the top to the bottom panels, respectively. 
The vertical lines indicate the change of Obs filter and the 
adopted Obs filter is labeled in the corresponding redshift 
range. 



magnitudes to derive k-corrections. We use a set of 
templates generated with the galaxy evolution model 
PEGASE.2 HFioc fc Rocca-Volmerange 1997| ). The tem- 
plates are computed for eight spectral classes including 
elliptical, spiral, irregular and starburst, with the initial 



mass function (IMF) from Rana & Basu (1992) with ages 
varying between 10 Myr and 14 Gyr. Dust extinction and 
metal effects are included, depending on the evolution sce- 
nario. We derive the absolute magnitude in the reference 
band Ref from the apparent magnitude in the band Obs: 



-.Obs 



DAI{z, Ho, ^ra, ^a) - KC{z, SED)(A.l) 



with DM the distance modulus and KC defined as follows: 



ifC(z,SED) = (fc^"^(z) + 



,Obs 



(^) 



(z))SED (A.2) 



with k^^-^ the k-correction in the reference band 
( |Oke fc Sandage 1968^ . To hmit the template dependency, 
the Obs band is chosen automatically to be the closest as 
possible to the Ref band redshifted in the observer frame. 
We do not correct the absolute magnitudes for the internal 
dust extinction related to the considered galaxies. 

To check the robustness of our absolute mag- 
nitude estimate, we use the GALICS simulations 
IIHatton et al. 2003|l . We extract a simulated catalogue 
with B, V, R, I apparent magnitudes and redshifts from 
the GALICS/MOMAF database. We apply exactly the 



16 



Ilbert O. et al.: Evolution of the galaxy global VVDS LF up to 2 = 



2 



method described before, to rederive the absolute magni- 
tudes. Fig lA. II shows the difference between our measure- 
ments and the 'true' absolute magnitudes from GALICS. 
When our procedure to limit the template dependency is 
efficient, the dispersion remains very small in comparison 
to the photometric errors. For instance, a ^ 0.03 — 0.04 
in the U-band, and we limit the template dependency up 
to z = 2. If our procedure can not be applied (since NIR 
data are not considered here), the dispersion increases. 
For instance, in the /-band the dispersion due to the k- 
correction is ~ 0.2. 

A. 2. Luminosity function estimators 

We describe in this subsection the four standard estima- 
tors implemented in ALF, the 1/Vmax, C+, SWML and 
STY estimators. 

A. 2.1. The 1/Vinax estimator 

The 1 /Vmax LF estimator (jSchmidt 1968|l is the most of- 
ten used to derive the LF, because of its simphcity. This 
estimator requires no assumption on the luminosity distri- 
bution. The 1 /Vmax gives directly the normalization of the 
LF, assuming implicitly an uniform spatial distribution of 
the galaxies. 

We consider a sample selected between bright and faint 
apparent magnitude limits, rrih and to/ respectively. The 
maximum observable comoving volume in which galaxy i 
can be detected is given by 

VobsA =11 - — —dujdz, (A. 3) 



where iv is the effective solid angle of the survey, and V is 
the comoving volume. Zmin,i and Zmax,i are the lower and 
upper redshift limits within a galaxy i can be included 
in the sample. The LF, (f>{M), is discretized in absolute 
magnitudes 



k=l 



where the window function W is defined as, 
1 ff -dM/2 <x < dM/2 



W{x) = 



otherwise. 



(A.4) 



(A.5) 



A. 2. 2. The C+ estimator 



|Lynden-Bell [ 1971) derived the non-parametric C~ 
method to overcome the assumption of a uniform galaxy 
distribution. We use a modified version of the C~, called 
0+ IjZucca et al. 1997ll . This method is based on the 
equality: 



d4 



dX 
C+'' 



(A.8 



with ip the cumulative luminosity function, d^p the vari- 
ation of the cumulative luminosity function between M 
and M + dM . dX is the number of observed galaxies be- 
tween M and M + dM and C+ is the number of galaxies 
brighter than M+dM with a redshift lower than the maxi- 
mum redshift observable. We use a sample sorted from the 
faintest to the brightest absolute magnitude. We note 
the value of C+ for a galaxy i. We introduce the weight 
Wi in Ct as follows: 



E 

j=l,Zj<Z„ 



(A.9) 



We note fi the variation of the cumulative luminosity 
function in the neighborhood of the galaxy i, between 
Mr and M+ . We can write the cumulative luminosity 
function: 

0=i j=l 



(A.IO) 



We impose the limit values /o = and V'o = 1 to nor- 
malize the cumulative luminosity function. We obtain the 
recurrence relation, used to derive the contribution of all 
galaxies in the sample: 



fi — ^+ 



ill 1 ^ S}=i fj 



The LF is given by rebinning the contributions J^: 



(A.ll) 



bkdM = Y f^W^W(Mk - Mi). 



(A.12) 



Poisson errors are associated as done for the 1 /Vmax esti- 
mator fEa. lA.7|l . 



is derived in each absolute magnitude bin k as follows: A. 2. 3. The STY and SWML estimators 



hdM = Y 



Vob 



-W{Mk - M,[ 



(A.6) 



where Ng is the number of galaxies of the sample and Wi 
is the weight applied to correct the unidentified sources in 
the field (see Section 3). We associate the Poisson errors 
to the 1/Vmax following IjMarshall 1985|l : 



The STY dSandage et al. 1979| ) and the SWML 
IjEfstathiou et al. 19^ hereafter EEP88) estimators 
are both derived from maximum likelihood method. The 
likelihood C is the joint probability of observing the 
galaxy sample, taking into account the observational 
selection effects. The principle of the SWML and STY is 
to maximize C with respect to the LF. C is given by: 



\ 



E 



N„ 



W{Mk - Mi) 



Vo 



bsi 



(A.7) C^l[p{M,,. 



n 



C/r m)dM j 



(t>(Mi) 



(A.13) 
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where Mf^i and Mf,^,; are the faint and bright observable 
absolute magnitudes of a galaxy i at redshift Zi. The 
weight is introduced in C following Zucca et al. (1994) 



The second derivative of the likelihood is given by: 



This weight artificially decreases the size of the er- 
ror contours derived from the analysis of C, then we 
balance this weight by the average weight < w >. The av- 
erage weight < w > does not affect the minimization of C. 

The STY assumes a functional form for the luminos- 
ity distribution. We use the empirical Schechter function 
HScher.hter 197611: 



S'^lnC ^ ]_\^^si wi5,jW(Mi-Mj)dM^ 

<w> 2^1=1 (0jdM/(g+l))2 
1 sr^Ng widM^Hl 



<w> 



= 1 (E"=r 't>kdMH2/(g+l)r 



(A.19) 



L* 



<j){L)dL = (t)*e-- 
The likehhood (ea \KJ^ can be written as: 



(A.14) 



with HI = H{Mb,i - Mi)H{Mr ~ Mf^i)H{Mb,i ~ 
M,)H{M,-Mf,i) &nAH2 = HiMi„i - Mk)H{Mk - Mf^i). 
The error bars of the LF (for a normalization given by the 
constraint) are given by the square root of the diagonal 
values of the covariance matrix. 

A. 2. 4. Luminosity function normalization 



InC 



{l/L*)j:fj,w,U-{a + l)lnL*j:l 



Na 



wMnc. + 1, ^) - T{a + 1, ^))] 



The estimators independent of the spatial density distri- 
bution (SWML, STY and C+) lose their normahzation 
while the normalization is directly done for the l/Vmax 
estimator. We adopt the EEP88 density estimator to re- 
(A.15)cover their normalization. The density n is simply the sum 
over all the galaxy sample of the inverse of the selection 
function: 



with r the incomplete Euler gamma function. We 
use the MINUIT package of the CERN library 
HJames fc Roo s 1995) to minimize -2lnC (MIGRAD 
procedure), to obtain the non-parabolic error for each 
parameter (MINOS procedure) and the error contour 
a - M* (MNCONT procedure). The crosses of the hke- 
lihood surface with InCmax — ^InC is used to compute 
the errors. The threshold AlnC is chosen in a standard 
way that depends on the desired confidence level in 
the estimate ( e.g., 2AlnC = 2.3 and 2AlnC = 4.61 to 
estimate the a — M* error contours with 68% and 90% 
confidence level ; 2AlnC = 1 to estimate the one sigma 
error for one parameter). 



1 



Vt 



total 



(A.20) 



The comparison with the l/Vmax normalization is a di- 
rect and independent check of the LF normalization. 
The parameter (/)* is directly related to the density 
0* (j}{M)dM = n. 0* is a function of a and M*. To 
estimate the error on (jj* , we derive (/)* for the extreme val- 
ues of the a — M* error contour at one sigma confidence 
level. We adopt Poisson errors when larger. 



The SWML does not assume any functional form for 
the luminosity distribution. The LF is discretized in ab- 
solute magnitude bins like the l/Vmax (see Eg. IA.4|) . We 
maximize InC with respect to 0fc to obtain the recurrence 
equation: 



with 



Ng WiH{Mb,i-Mj)H{Mj-Mf^i) 
~^ J2k=T <t>hdMH(Mh^i-Mu)H{Mk-M}^i) 



(A.16) 



if a; < -dM/2 
H{x) = { x/dM + 1/2 if -dM/2 <x< dM/2 

1 otherwise. 



(A.17) 



We add a constraint g on cfik and rewrite the likelihood 
as InC = InC -\- Xg where A is a Lagrangian multiplier. 
Following EEP88, we choose g = J2k=i <l>kdM {Lk / L jf - 
1 with Lf the fiducial luminosity and (3 a constant. The 
error bars are derived from the covariance matrix, denoted 
C, defined as the inverse of the information matrix /: 



